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JiJ We point out that there is a nonabehan instabihty for a nonabehan plasma which does not 

^^ allow both for a net nonzero color charge and the existence of field configurations which 

^ — I 
;;• are coherent over a volume v whose size is determined by the chemical potential. The basic 

\^ process which leads to this result is the Schwinger decay of chromoelectric fields, for the 

;_ I case where the field arises from commutators of constant potentials, rather than as the 

curl of spacetime dependent potentials. In terms of the fields, instability is obtained when 

Tr(D°F^^ D^F^,) > 0. 



The identification of the deconfined phase of quarks and gluons at the Relativistic Heavy 
Ion Colhder, a phase akin to a nonabelian plasma, has led to a number of investigations 
on instabilities in a nonabelian plasma [T, "2]. While some of these are concerned about 
an upgraded version of instabilities in an abelian plasma, such as the Weibel instability, 
there have been numerical studies of the evolution of instabilities in the hard thermal loop 
approximation and beyond. The purpose of this note is to point out that there is an 
instability, and a certain no-go statement, which is quite general and arises purely from 
nonabelian effects. It is fairly straightforward to understand how this effect arises. For a 
statistical distribution of nonzero color charge, we need a chemical potential. Because the 
charge is nonabelian in nature, the chemical potential is a matrix in the Lie algebra of the 
color group. In fact, it may be viewed as a background value for the time-component of the 
potential ^o = —if^Afj, where t" form an orthonormal basis for the Lie algebra of the color 
group G. (We may actually take this matrix to be diagonal, but it is not important at this 
stage.) If we have a constant background Aq, then there is an electric field generated via 
the commutator term [Aq, Aj] in the field strength tensor. For modes of Ai of wavelength 
A, this gives an electric field approximately constant over this length scale. This electric 
field will then develop a Schwinger instability decaying via pair production. If the particles 
which are produced have a mass, there is an exponential suppression, but in the nonabelian 
plasma, we have effectively massless modes. The end result of this argument is the following. 
Consider the plasma coarse-grained over a distance scale A. Then one possibility is that 
the color charge density is zero when coarse-grained over this scale. The other possibility 
is that the plasma cannot have Ai which are coherent over length scales exceeding A. This 
is the essence of our no-go statement. 

The possibility of color charge density being zero has been studied in the context of 
color superconductivity [3J . In the limit of large baryon number density, we expect a color 
superconducting phase and it is important to have color neutrality. Such a requirement 
can be imposed on analyses of color superconductivity, but how it is achieved is really 
a dynamical issue. (This is not the setting for our question. We are concerned about 
a deconfined state, not superconducting and for us the baryon chemical potential can be 
zero. But there are points of connection.) Nonzero charges can lead to large electric fields 
which are unstable, can lead to energy being nonextensive and this is one reason why stable 
matter must be neutral under gauge charges [3]. Nevertheless, it is interesting to analyze 
some of the nuances of how neutrality is achieved. Since the chemical potential may be 
taken as a background value for Aq, the corresponding equation of motion (or integration 
of the constant mode of Aq in the functional integral) seems to imply zero color charge. 
Strictly speaking this argument needs to be qualified, since it is equivalent to imposing 



the Gauss law integrated over functions which do not vanish at spatial infinity. The true 
gauge transformations of the theory go to the identity element at spatial infinity and so test 
functions for the Gauss law must vanish at infinity. Imposing the Gauss law with constant 
values for the gauge parameters is equivalent to eliminating all charged states by fiat, which 
we do not want to do. One can use a compact spatial manifold and then approach the limit 
of large volumes to preserve the zero charge condition. This provides a method for carrying 
out the analyses, including many of the calculations in the literature, but it is not quite 
an explanation. All this makes it useful to ask the question we are asking: If we have a 
deconfined state of gluons (and may be quarks), and we try to have nonzero color charge, 
what instabilities can arise? 

The density matrix for a statistical distribution in equilibrium is given hy p = exp[—{H— 
^i fJ-i Qi)/T\ where H is the Hamiltonian, Qi are conserved charges, pi are the correspond- 
ing chemical potentials and T is the temperature. We are interested in time-dependent 
processes in this distribution, so we are concerned with real-time propagators and vertices 
averaged over states with the density matrix p. The result is equivalent to calculations at 
zero chemical potential, but with a Hamiltonian H — ^^PiQi- Since the constant mode 
of Aq couples to (5,it is clear that we can treat /x as a background value for Aq. Consider 
now the nonabelian charge density due to quarks, say, Jq = qj^f^q, or its matrix ver- 
sion, (Jo)ij = Qi^^qj, i,j being color labels for the quarks. Under a gauge transformation 
g{x) G G, this matrix changes as 

^0 ^ J^ = g'^ Jo 9 (1) 

It is thus possible to choose g{x) such that Jq is diagonal at each point. In other words, the 
gauge-invariant information contained in Jq may be taken as the diagonal charge densities. 
Thus, to specify a charge distribution, we need only chemical potentials for the Cartan 
elements of the Lie algebra. There are other ways to see this as well. For example, if 
the charged particles form some irreducible representation R (which may be thought of 
as arising from the decomposition of a product of the representations of the individual 
particles), then we know that such a representation can be obtained by quantizing the 
co-adjoint orbit action 

S = i fdtY,WkTr{hkg-^g) (2) 

k 

where Wk are the highest weights defining the representation R and hk are the diagonal 
generators of the Lie algebra. We see that the diagonal charges are sufficient for our purpose. 
In a statistical distribution, we have to think of such a representation for the global color 
charge over each coarse-grained volume element, and this action can be generalized to obtain 
the fluid flow equations for color charge [IJ. 



In the case of a nonabelian plasma, there is an added comphcation. While it is possible 
to define a gauge-covariant charge density for the quarks (and other matter particles), there 
is no gauge-covariant charge density for the gluons. The integrated total charge has a 
gauge-invariant expression. The chemical potential, introduced as a background value for 
^0 does couple to this global charge correctly. This also leads to terms quadratic in fi in 
the action, which is to be expected since the current for a charged bosonic system depends 
on A^ in addition to the charged fields themselves. All these effects are included in the 
replacement Aq ^ Aq + fi. Since the diagonalization of the charge density happens only by 
choice of g{x), the general ansatz for the background value of Aq is 

Ao = g'-^fi g + g'^dog (3) 

The group element g can be removed by an overall gauge transformation, 

^0 -^ g^o9'^ - dogg^^ = /x 

Ai -^ gAig-^ - digg-^ (4) 

Designating the new spatial components of the potential as Ai again, we see that we can 
use ^ as the background value for Aq. 

Calculating the effective Lagrangian 

We shall carry out the calculations in Euclidean space. While this is not necessary, as 
for many other calculations at finite temperature and density, this is slightly simpler. This 
means that the background value of the Aq becomes imaginary. Thus the basic calculation 
to check for instability reduces to the following. Taking constant matrices for Aq and Ai as 
the background values, we consider fluctuations in the fields. The integration of the action 
to quadratic order in the fluctuations leads to the standard determinant. This has to be 
evaluated as a function of the background values. The result is then analytically continued 
to imaginary values of the background ^o- The result can then be analyzed for instbilities. 
The instability of interest to us is the Schwinger decay of the chromoelectric fleld. This 
has been studied in some detail in the nonabelian case for electric flelds which are given by 
the curl of the gauge potentials [5j, but, here, we are interested in the case when the fleld 
arises from the commutator term of the potentials. For the calculations which follow, we 
will consider the group SU{2) since it is sufficient to capture the effect we are interested in. 

The integration over the quadratic fluctuations can be phrased as an effective Lagrangian 
given by 
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where the second term is the contribution from the ghosts. Here D^ = (5^ + Ap){d^ + A^^) 
is the gauge-covariant Laplacian with the background field A't; it is a 3 x 3-matrix in color 
space, as indicated by the color indices a, 6. Thus the operator —{D'^Y^r]^^ — 2/°'^''F^^ 
can be considered as a 12 x 12-matrix, in addition to its coordinate space properties. The 
evaluation of the action will follow a method which is similar to what was used many years 
ago by Brown and Weisberger [6|. Writing the SU{2) field Afj^' = f°-'^^A'l 
simplify D^ as 






- {D')ab =p' + Y~Yab~ 2ip ■ A,b (6) 

where p^ = —idfj_, Y°-^ = A°^A^ and Y = TrY"-^. The matrix Y"-'' can be diagonalized by 
a suitable gauge transformation, with eigenvalues A^. These eigenvalues give the gauge- 
invariant characterization of the chromoelectric and chromomagnetic fields. The A's are 
positive in the case of Euclidean signature for the contraction of spacetime indices in A"j^A , 
but one eigenvalue can be negative with Minkowski signature. In the Euclidean metric we 
are using, we can always choose 
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Aa5^, 



a,n = 1, 2,3 



(7) 



With this choice 



Yab + 2iip- A)ab 



Ai -2^p3^/}^ 2ip2\^ 

'^ipsVh A2 -2^pl^/X^ 
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(8) 



For our purpose, it is not necessary to consider this matrix in full generality, we can take 
A3 = 0. In this case the only nontrivial component of the field strength tensor is F12 = 
—F21 = \/AiA2. In this case, schematically, we have 



[Yab + 2i(p • A)ab] r]^u + 2Fabf_Lu 



Y + 2ip-A 2Fi2 

-2Fi2 Y + 2ip- A 

Y + 2ip-A 



















Y + 2ip- A 



(9) 



where each block is a 3 x 3 matrix in color space. From this block diagonal form, 



Tri2xi2 exp [s{{Y + 2ip ■ A)r]^u + 2F^^}] 



2 Tr3x3 e^(^+2^P-^) 

+ Trexe e^[(^+2*P-^)'?.^+2^."] (10) 



The first term on the right hand side cancels exactly the similar contribution from ghosts. 
The remaining 6x6 matrix corresponds to the indices 1,2, for spacetime and the 3x3 
matrix in color space. The effective Lagrangian is thus 

1 f d^p f°° ds 
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For evaluating the remaining trace, it is convenient to use the integral representation 
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^l^e— ^logdet(z-X) 
2tti oz 



(13) 



where the integration contour encircles all zeros of det(z — X). 
The determinant is easy to evaluate, 
det(z -X) = {z^ + z\Xi + A2) - ['ipHzXi + A?) + 4pl{zX2 + A^) + 3ZA1A2] }' 



Iz[z^ + z{li + l2)-3lil2 



Apl{zh + ll) + Apl{zh + ll) 



(14a) 

Z[z^ + z{h+l2)-2,lil2] Ji ^-^^^^ 

When this is used in ( 11|13 ), with the dz carried out, we get contributions from the poles 
which correspond to the roots of the cubic polynomial inside the braces in (14a). It is then 
convenient to split the expression for ie// as Li + L2 with 

d'^p /■^ ds 
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The evaluation of Li is simple. The zeros of the relevant cubic polynomial are z = 
and z = z± with 



1 



z± 
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We then find 
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r is the Eulerian gamma function. Notice that there are singularities in this expression for 
D = A. These are, of course, the standard renormalization singularities and can be isolated 
by expanding (/i)^~^Li in powers of e with D = A — e. (The ^-factor is the usual one for 
ensuring the correct dimension for Li.) This leads to the expression 
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where fi? = 47re~'^/x2, 7 being the Euler-Mascheroni constant. 



+ Oie) 



(18) 



The first term on the right hand side of (18) is the potentially divergent contribution 



which is removed by renormalization. The remainder gives the finite expression we need for 
Li. 

The evaluation of L2 is a little more involved and is sketched out in the appendix. The 
final result is 
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(19) 



(20) 



and A2 is obtained by the exchange /i -H- I2 in the above expression. In (19) also, there is 



a potentially divergent contribution arising from the large z behavior of the integrand. Its 



removal, along with the potentially divergent term from ( 18 ) is discussed in the appendix. 

The nature of the instability 

We are now in a position to consider how instabilities can arise from these results. In 
continuing the expressions for Li, L2 to Minkowski space, one of the directions has to be 



identified as the time-direction. We will take this to be the 1-direction. Chromoelectric 
fields in Minkowski space will thus correspond to the choice /i < 0, I2 > 0. The choice 
of /i,/2 > will correspond to the purely chroniomagnetic case, with 1-direction being 
interpreted as spatial direction now. We will consider various possibilities for the A's one 
by one. 

Case a : 

Consider first the case of li < 0, I2 > 0, h + h > 0. In this case, the factor {li + ^2)^ + I2I1I2 

is positive for I2 ^ |/i|. For this region 

,, + ;. + ., = 'i±ii±vlTiF±Mi >„ pi) 

and hence there is no instability in Li. As we come down in the value of I2, this factor 
changes sign at h = {7 + \/48)|/i|. For the region |/i| < /2 < (7 + \/48)|/i|, the quantities 
h + h + -2+ and h + I2 + Z- are complex conjugates of each other. Writing these as ae^*^. 



we can easily see from (18) that there is no imaginary part in L\ for this region as well. 



Thus, there is no instability resulting from L\. 

Turning to ImL2, notice that we can set Z) = 4 at this stage because the integration 
range for z for the imaginary part does not extend to infinity and so the issue of divergences 
do not arise. The analysis of L2 then reduces to the analysis of the condition A2{^z) > 1. 
The polynomial factor in the denominator of the A's, namely, that (z + l\){z + I2) — ^l\l2 = 
z^ + z{li + I2) + 3|/i|/2 is easily seen to be positive. Thus Ai{z) < and the factor \/l — xAi 
is real for the full range {z > 0) of integration for z. On other hand, ^42(2;), whose numerator 
is 4:zl2{z + li) will show a change of sign for z = —li > 0. However, even though A2{z) > 
for z > li, we have A2{z) < 1. This is easily seen from the fact that 

{Z + I1 + h) [{z + h){z + I2) - 4/1/2] >{z- \h\) [{z - \h\){z + h) + 4|/i| h] (22) 

The quantity in the square brackets on the right hand side is > Azl2 for z > \li\. Thus the 
factor \/Y^^~xA2 is also real and hence there is no instability for this case from either Li or 

L2. 

Case b : 

Now we turn to the case h < 0, I2 > 0, h + 12 < 0. The region (7 — -v/48)|/i| < h < \li\ has 
complex conjugate values for I1+I2 + z± and there is no imaginary part resulting from the 
last two terms in Li, as in the previous case for I2 < {7 + ^48)1^1 1- There is an imaginary 
part from the log(/i + I2) term in Li, which will give an instability for this range of h- For 
l2<{7- \/48)|/i| (or |/i| > (7 + \/48)/2) we have 

li + I2 + z± = < (23) 



200^ 



15 20 



Figure 1: Sample graphs of Ai{z) for 1 < |/i|//2 < 7 + \/48 (left) and for \h\/l2> 7 + \/48 
(right). The value of Ai between 15 and 20 is large and positive and outside the frame of 
the graph on the right. 
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Figure 2: Sample graphs of ^2(2) for 1 < |/i|//2 < 7 + V48 (left) and for \h\/l2> 7 + V48 
(right). 



There is then a nontrivial imaginary part in Li which leads to an instability. Thus we get 
instability from Li for all h, h corresponding to this case. 

Turning to L2, we may notice that the factor (z + li + I2) in the denominator of Ai, 
A2 changes sign at z = —(h + h)- The additional factor in the denominator, namely, 
[{z + h){z + I2) — 4/1/2] has two positive roots if I/1I//2 > 7 + \/48 w 14. Otherwise, there 
are no real roots and this factor is positive. The graphs of Ai{z) as a function of z are as 
shown in Figfl] We see that for all values of I/1I//2, there are regions of z-integration for 
which Ai{z) > 1, leading to an imaginary part for L2. Similar statements apply for A2, see 
FigjSj 



Case c : 

Even though it is not germane to our present discussion, we may note that if we have the 

purely chromomagnetic case with /i > 0, ^2 > 0, then 

/i + ^2 + ^- = J — < (24) 



Thus the last term on the right hand side in (18) has an imaginary component. For L2, the 



polynomial (z + li)(z + l2) —4/1/2 in the denominators of ^1, A2 has roots z±. For li,l2 > 0, 
one root is negative and the other is positive. Ai{z) is positive for z > z-|_ and goes to 
zero for large z, with Ai{z) — t- 00 for z — Z-^- — t- 0+. Thus there is a range of z for which 
\/l — xAi has an imaginary part. Again a similar statement applies to A2- Thus for both 
Li and L2 we get an instability for li,l2 > 0. This is the well-known vacuum instability in 
a chromomagnetic field. 

It is interesting to characterize the instability in terms of invariants of the field. We see 
easily that F^^^F"'''' = {TrYf - TrY^ = 2/1/2, {D''F^"')''{DaFf,^y = 2/i/2(/i + h). We may 
then summarize our results as 

. „ „„ , I > Instability , , 

[ < No instability 

Comments 

The instability we are discussing is quite general and hints at how statistical distributions 
tend to move to color neutrality or a disordered state with no coherent fields over distances 
long compared to the dimension given by the chemical potential. The calculation itself may 
be taken as the derivation of Schwinger decay of chromoelectric fields for the case when the 
field is generated by the commutator term, rather than the curl of the potentials. 

This work was supported by U.S. National Science Foundation grant PHY-0855515 and 
by a PSC-CUNY award. 

APPENDIX 

Calculation of L2 

For L2, we start with the representation 

logA= %-t-e-'^) (26) 

Jo T^ 
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Figure 3: Deformation of contour 
for branch cut at z = —t — li — I2 



Figure 4: Contour for evaluating L2 



Using this and eUminating dz by partial integration, the expression (15b) for L2 becomes 
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(27) 



(28) 



for k = 1,2. For the second line of equation (27) we have carried out the p-integration. 



Note that the exponents involving pj, show that we need o take the z-contour to be large 
enough, \z\ > 2\/tl. Effectively, this means that we should do the z-interal before doing the 



(29) 



i- integral. In (|27|), we can further carry out the s-integration to get 
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The factor (z + f + /i + 12)^^^^'"^ shows that, for the z-integration, we have a branch cut 
along the negative real axis starting at z = —t — li — l2- We can deform the original contour 
which is a large circle around the origin, via the contour shown in FigjSj to the contour in 
Fig|4j Notice that because of the arguments given earlier, the branch point z = —t — li — l2 
is always outside the original contour, while the singularities of the square root factors are 
always inside the contour. Integration along the cut in Figj4] gives 
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where 
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(^ + /l + /2)[(^ + /l)(2 + /2)- 4/1/2] 

and A2 is given by the same expression with li 0/2- Changing the order of integration and 
making the substitution t = zx, we finally get 



00 /"i 



L2 = —-, — tfttt; — -, r^ / dz I dx 

(4^)^/2r(Z)/2) Jo Jo 



(1-^y 



'1+15/2 



-l+D/2 



X 



1 



1 



^/l - xAi^/1 - xA2 



(34) 



This is the expression quoted in the text. 
Renormalization : A consistency check 



The potentially divergent part of Li was obtained in equation ( 18 ) as 
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Using the expressions for z± from (16), this simplifies to 
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The expression for L2 can be recast as 
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where r = l/z and 
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and A's correspond to A's with 2; = 1/r; i.e., 

4T/i(l+r/2) 
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[1 + r(/i + /2)] [(1 + t/i)(1 + T/2) - 4t2/i/2] 
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(39) 



with li <-)• I2 to obtain A2 from Ai. The divergence now corresponds to small values of r. 
Carrying out a small r-expansion, 



G(r) 
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[Z)/i/2 + p-l) (/? + /!)] + O(t') (40) 



We can use this expansion in (37) and integrate; we are interested in small r region, so we 



use a cutoff e""^ in the integrand. (Whether we use this or something else, such as e~"'^ for 
some a does not matter for the term of the form r((4 — D)/2).) The term proportional to 
1/e is then found to be 
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Combining this with ( 36 ) , we find 
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This is the expected and correct renormalization of the action, and is consistent with the 
/3-function of 

for SU{2). 
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